Quantum to Classical Transition from the Cosmic Background Radiation 



M.C. de 01iveira\ N.G. Almeida^, S.S. Mizrahi^, and M.H.Y. Moussa^ 
^Department of Physics, University of Queensland, QLD 4 072, Brisbane, Australia. 
^ Departamento de Fisica, Universidade Federal de Sao Carlos 
Via Washington Luis, km 235, Sao Carlos 13565-905, SP, Brasil 
(February 1, 2008) 

We have revisited the Ghirardi-Rimini- Weber- Pearle (GRWP) approach for continuous dynamical 
evolution of the state vector for a macroscopic object. Our main concern is to recover the decoupling 
of the state vector dynamics for the center-of-mass (CM) and internal motion, as in the GRWP 
model, but within the framework of the standard cosmology. In this connection we have taken 
the opposite direction of the GRWP argument, that the cosmic background radiation (CBR) has 
originated from a fundamental stochastic hitting process. We assume the CBR as a clue of the 
Big Bang, playing a main role in the decoupling of the state vector dynamics of the CM and 
internal motion. In our model, instead of describing a continuous spontaneous localization (CSL) 
of a system of massive particles as proposed by Ghirardi, Pearle and Rimini, the Ito stochastic 
equation accounts for the intervention of the CBR on the system of particles. Essentially, this 
approach leads to a pre-master equation for both the CBR and particles degrees of freedom. The 
violation of the principle of energy conservation characteristic of the CSL model is avoided as well 
as the additional assumption on the size of the GRWP's localization width necessary to reach 
the decoupling between the collective and internal motions. Moreover, realistic estimation for the 
decoherence time, exhibiting an interesting dependence on the CBR temperature, is obtained. From 
the formula for the decoherence time it is possible to analyze the transition from micro to macro 
dynamics in both the early hot Universe and the nowadays cold one. The entropy of the system 
under decoherence is analyzed and the emergent 'pointer basis' is discussed. In spite of not having 
imposed a privileged basis, in our model the position still emerges as the preferred observable as in 
the CSL model. 

PACS number(s); 03.65.Bz, 32.80.-t, 42.50.Dv 



I. INTRODUCTION 

In the last decade several proposals to modify the standard Hamiltonian dynamics, ranging from master equations 
to stochastic quantum mechanics, have been advanced to try to set up an unified description for microscopic and 
macroscopic physical phenomena. In the pioneer work by Ghirardi, Rimini, and Weber quantum mechanics with 
spontaneous localization (QMSL), the state vector collapse, leading from quantum to classical dynamics results from the 
instantaneous action of a spontaneous random hitting process. Such a Poisson process is described by a "localization" 
operator, a gaussian function acting on each microscopic constituent of any system. The localization operator carries 
two free parameters; a mean frequency A and a localization width a"^/^, understood as new constants of nature (the 
spontaneous localization is argued to be a fundamental physical process). Through these basic assumptions the QMSL 
consists in an explicit model allowing an unified description for microscopic and macroscopic systems. It forbids the 
occurrence of linear superposition of states localized in far away spatial regions and induces a dynamics that agree 
with the predictions of classical mechanics. 

Pursuing the program of the QMSL model, Diosi presented an interesting connection between the original 
GRW hitting process and a modified Schrodinger equation. Another significant achievement concerning a dynamical 
reduction model, a stochastic equation for physical ensemble, was reported by Gisin Next, Pearle [Q described 
the QMSL model through an Ito stochastic differential equation. Basically, Pearle replaced the Poisson process of 
instantaneous hits in GRW model by a Markov process described as a stochastic modification of the Schrodinger 
equation, so that a continuous evolution of the state vector was accomplished. By considering a specific choice of 
the operators defining the Markov process (expressed in terms of creation and annihilation operators), Ghirardi, 
Pearle, and Rimini [|| have described the mechanism known as continuous spontaneous localization (CSL) of systems 
of identical particles (the QMSL model has consistency only in the case of systems of distinguishable particles). 

Other investigations dealing with dynamical reduction models have recently been considered among them 
it is worth to mention the model for intrinsic decoherence proposed by Milburn |^]. While in the GRWP model 
the addition of stochastic terms in the Schrodinger evolution automatically destroys the quantum coherence of the 
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physical properties of the system that attain a macroscopic level, the modification of the Liouville equation proposed 
by Milburn destroys the coherence even at microscopic level. 

In the CSL model the Ito stochastic equation for the evolution of the state vector reads 



d\^) = ^—f:Hdt + dh- ^idh)^j (1) 

where dh is a linear self-adjoint operator, whose random fluctuation may increase or decrease the norm of the state 
vector. Using Ito formula (with the notation \dip) = d\tp)), 



d II V f = ii^m + {di^w + wi#), (2) 

it is easy to see that Eq. (0) does not conserve the norm of 1-0). Thus, the introduction of a norm conserving nonlinear 
process is mandatory. This process, whose random operator depends on the state vector, reads 

d\cl>) = (^-^Hdt + dh^ - ^Jdh^^ 10). (3) 

Now, it is necessary to distinguish between raw ( Eq. (|l|)) and physical ( Eq. (^)) ensembles of state vectors to correctly 
understand the effect of the non-Hamiltonian terms. To this end a precept is adopted, namely, that the square norm of 
each (unnormalized) state vector represents the weight associated with that (normalized) state vector in the ensemble 
coming from the Ito stochastic equation |^,||. This precept is a generalization of the GRW assumption that the 
frequency of hits is proportional to the squared norm of the state vector. Therefore, in the GRW prescription the 
quantum theory prediction about the associated probabilities in a measurement process is recovered. By considering 
such a precept for the physical ensemble, the linearity of the raw equation and the Markov nature of the Ito stochastic 
process leads to the physical stochastic differential equation for the A''-particle state vector 

d|*Ar) = (y-^Hdt + Z.dB - ]^-iZKZdt^ \^n), (4) 

where Z = {Zi} are operators on the Hilbert space of the system and the set of random operators B = {Bi\ is 
characterized through a real Wiener process, satisfying the following means and correlations over ensemble 



dB, 0, dBdBj = 7*ydt. (5) 



The statistical operator pn = \^n){^n\ of the physical ensemble and its evolution equation are directly obtained 
from Eq. (||); using the Ito calculus in evaluating dp^/dt one gets 

^ = [if, pn] + TZpAT.Zt - ^ {Zt.z, p^} , (6) 

which is exactly the Lindblad 1^ form for the generator of a quantum dynamical semigroup. 

In the present work our main concern is to achieve the decoupling between the state vector dynamics of the center- 
of-mass ( CM) and internal motion of a system of particles. In the GRWP model this decoupling results from a 
hypothesis of spontaneous localization of the system's wave function due to a fundamental stochastic hitting process 
on the particles, which induces an increase of total mean energy of the Universe claimed to be the origin of the 
Cosmic Background Radiation (GBR). Gontrarily to this argument, in the present work we assume the point of view 
of standard cosmology: the nowadays GBR is a clue that the Universe began its expansion from a Big Bang [ p^ . 
This assumption is introduced with the purpose to avoid the unconventional increase of the total mean energy of 
the Universe. Formally, we hypothesize that the state vector, the Hamiltonian H and operators Z, Z^ in Eq.(U) 
represent both, the system of particle and GBR radiation; the set of random functions {Bi\ describes the intervention 
of the GBR on the system and substitute the spontaneous localization process. Instead of elaborating on the formal 
microscopic problem of the interaction of a system with a reservoir |^, we assume ad-hoc that the evolution of 
the system of particles, under the influence of the GBR, is described by an Ito equation having stochastic coupling 
parameters. 

Therefore, in the present conservative continuous reduction model (the total energy of system plus GBR is conserved) 
we argue that: 1) the increase or decrease of the system's mean energy is attributed to the GBR; 2) the positional 
space is not privileged with respect to the momentum space, as required when the localization operator is involved; 
3) we do not claim for an additional assumption to decouple the collective and internal motion, namely the width 
parameter a^^/^ ~ lO^^cm in the GSL model; 4) as above-mentioned, more admissible results are obtained for 
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decoherence times, while in the CSL model the value 10 s obtained for a system of particles to undergo from 
quantum to classical dynamics seems to be too large (as well as the localization width a~^^'^ ~ 10~^cto also seems too 
large when considering typical atomic distances about 10~^cm, or even superposition of the center-of-mass coordinate 
different by more than about a^^^^ ^^), and finally, 5) instead of the two free parameters required in the GRWP 
model (a^^/^ and the mean frequency A), the random function describing the interaction between the system and 
the CBR carries just a single strength parameter with dimension of inverse of time. In fact, the coupling constant 
of the CBR photons to the TV-particle system, as the strength parameter in the GRWP model, defines the inverse of 
a characteristic time which is associated to the net effect of the random pseudo- "potential" dh Also, as in the 
GRWP model, our strength parameter is such small that nothing changes in the Hamiltonian dynamics of a single 
particle even in the case in which it has an extended wave function 

Finally, we mention that Joos and Zch , have previously argued that scattering of photons even at a relatively 
low temperature can induce the localization of the wave packet of a macroscopic system. So, their treatment, based 
on a master equation proposed by Wigner |Q , suggests that the intergalactic cold CBR cannot simply be neglected 
[ psf . The model here presented goes exactly on this point, i.e., we consider the process of random scattering of the 
CBR photons by a system of particles as responsible for the superselection rules and the micro to macro transition of 
its dynamical description. In this way, despite inducing the superselection rules the CBR also induces the mechanism 
of separating the center-of-mass (CM) coordinate from the internal motion. Besides, we present a brief cosmological 
analysis of our results, discussing the roles played by both the CBR temperature and the number of particles of the 
system, in its way from quantum to classical dynamics, as the universe evolved from a hot to a cold state. 

In Section II we briefly review the GRWP model presenting its main achievements. In Section III we construct our 
model: beginning from an Ito stochastic equation we derive a pre-master equation for a system of N particles and 
the CBR; tracing over the CBR degrees of freedom we obtain a master equation for the system of particles only and 
In Section IV we show that structurally it shows exactly the Lindblad form. In Section V we estimate the coupling 
parameter and in Section VI we estimate the decoherence time for the system of particles. In Section VII we show 
that at low temperature limit our master equation and the GRWP Ito equation are equivalent, thus this last one is a 
particular situation of the former; these equations allow the decoupling of the state vector dynamics into two separate 
equations, one for the CM and the other for the internal motion. In section VIII we calculate the entropy and analyze 
the problem of selection of a preferred basis. Finally, in Section IX we present a summary and conclusions. 



II. THE GHIRARDI-RIMINI-WEBER-PEARLE MODEL 



As explained in the introduction, in CSL model the random operator dh contains in its definition the length 
parameter a~^^^ and a strength parameter ^ which is related to the mean hitting frequency A. In this section we 
present a brief review of the CSL model as a class of Markov processes in Hilbert space [gj. We will consider a 
system of N identical particles so that the localization operator must involve globally the whole set of particles in 
order to preserve the symmetry properties of the wave function |ll6|| . For this purpose let us consider the creation and 
annihilation field operators a^(q, s), a(q, s) of a particle at the point q in some reference frame with spin component s, 
satisfying the canonical commutation or anticommutation relations. From these operators a locally averaged number 
density operator is defined as 



iV(x) = (£)3/2^ J d^qexp 



--a(q-x)2 



at(q,s)a(q,s). 



(7) 



The operator N{yi) is self adjoint and its commutator for different values of x vanishes. The total number operator 
is defined a,s N = J d^xA^(x), and the symmetrized (antisymmetrized) states containing n particles at the indicated 
positions, 

I q, s) = A/'at(qi, si)...a\q^, s„) | 0), (8) 
constitutes the normalized common eigenvectors related to the eigenvalue equation N{x) \ s) — n^x. \ q, s), with 



N 



--a(x- q,) 



(9) 



Applying |^,|6[ the stochastic process established by Eq. (^ to a system of identical particles and considering the 
locally averaged density operator defined by Eq. (|^), one gets the physical stochastic nonlinear differential equation 
for the state vector as 
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N 



-iHdt + J d^xN{x)dB{x) J d^xN^{yi)dt 



^n). (10) 



where the Wiener process i?(x) satisfies 



dB(x) = 0, (11a) 
dB{-K)dB{y) = C(5^(x - y)dt. (lib) 

From Eq. ( [To| ) the evolution equation of the iV-particlc statistical operator obtained from Ito calculus reads 

^ = -i[H, pn]+cJ d3x7V(x)p^iV(x) - 11 dWiK), . (12) 

and it can be checked that taking A = C('^/47r)'^/^, Eq. (|l2|) reduces to the correspondent equation for a single particle 
considered in the QMSL model. 

To discuss the physical implications of the modified dynamical equation (^0|) , the separation of the CM motion will 
be made. If Q is the CM coordinate of the system and its internal coordinates (measured from the CM of the 
particles), one can define the particle coordinates as 

q. -Q + q^({ri}), (13) 

where {ri} represents a set of 3A^ — 3 independent variables. In the GRWP model the set {ri} does not contain 
macroscopic variables. As a consequence, assuming that the Hamiltonian can be written as H — Hq + Hr^, we 
consider the wave function 

0(q,s) = *(Q)x(ri,s), (14a) 
x{ri,s)=(^^^A{rus), (14b) 

where the symbol (^) specify the symmetrization or antisymmetrization of the internal coordinate wave function. 
Under the assumption that the length parameter a-i/2 is such that the internal wave function A(ri,s) is sharply 
peaked around the value riQ of r (with respect to a"^/^ ), the action of the operator A^(x) on the wave function (14a) 
turns out to be 

iV(x)*(Q)x(ri, s) = F(Q - x)*(Q)x(ri, s), (15) 

with 



F(Q - x) = {^^" [Q + q«(^o) - ^]'} 



(16) 



Therefore, the operator -/V(x) acts only on ^> and the separately normalized wave functions ^ and x satisfy the 
equations 



iHqdt + J d^xF(Q - x)dB(x) J d^xF2(Q - x)dt 



^, (17a) 



dx = -iHr^xdt. (17b) 

By assuming a large enough length parameter and an internal wave function which is independent of the macroscopic 
variables, the internal motion decouples as in the absence of the stochastic terms in Eq. (^0|). From this fact, the 
reduction rates which are characteristic of the GRWP theory together with the position and momentum spreading 
can be obtained. In particular, in the positional representation of Eq. (|l2|), it is possible to verify with the help of 
the macroscopic density approximation and the sharp scanning approximation that the macroscopic frequency 
associated to the system of identical particles is 

T = CDonout- (18) 

Here a homogeneous macroscopic body of density Dq was considered and Uout is the number of particles of the body 
at position Q' that do not lie in the volume occupied by the body at position Q". In the case of distinguishable 
particles, one gets the direct result 
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Ac 



nA, 



(19) 



71 being the total number of particles, so that for a typical macroscopic number n « 10 , one obtains Xcm ~ 10 s, 
as mentioned above. 

The position and momentum spreading obtained from the approximations leading to Eq. (Esl), are written as 



where the suffix s indicates the Schrodinger evolution, and 



<5. 



(dF{y)-V 
V dy, ) 



(20a) 
(20b) 

(21) 



Now, using the macroscopic density approximation applied to the identical particles system, Eq. (Hq) is modified 



to 



F(Q-x)= j d3yD(y)(^)3/2exp 



--a(Q + y-x)2 



(22) 



where D{y) is the number of particles per unit volume in the neighborhood of the point y = Q + y. The evaluation 
of the factor 5i for the case of a homogeneous macroscopic box containing the N particles, through the Eq. (|2^) gives 
the result D 



5, = {a/irf'^DlS,, 



(23) 



where Si is t he tr ansversal section of the macroscopic box. 

From Eq. (20b) it is evident that the momentum variance implies that the CM energy increases per unit time as 



t 



M 



10"^^(5 cm s~^)S, cmr'^, 



(24) 



with the GRWP choice ^ iQ-s cm together with Dq ~ 10^^ cm ^ . From the requirement that the macroscopic 

frequency associated to the system of identical particles Eq. ( [l8| ) is exactly the same as for distinguishable particles 
Eq. (|l|), GRWP have chosen ( ^ 10"^° cm^s'^. 



III. DECOHERENCE FROM THE COSMIC BACKGROUND RADIATION 

Our approach uses the stochastic dynamical equation (|^), where we identify the continuous component (in frequency 
space) of the operator responsible for the interaction of the iV-particle system to the GBR as 



N 



Z(17) = ^(A(r!)at +At(f7)a^ 
fc=i 



(25) 



where 



and aj. is its hermitian conjugate ^ 



afc 



Ofe.i, Ofc/ 



1 



\j2hmijj 



(26) 



jjj, c — x,y, z j , qfc and pfc are respectively position and mo- 
mentum operators of the k*^ particle of mass m. ?ia; is a characteristic energy of the system of particles associated to 
the quantum fluctuation of the GM. The operators A^(yi), A{Vl) stand for the creation and annihilation of a quantum 
of energy Ml from the environment. The coupling parameter is defined by the continuous stochastic Wiener process 
B(ri) satisfying 
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dB{n) = 0, (27a) 
dB,{VL)dBj{VL') = -f{VL)5^^j5{Q. - n')dt, (27b) 

with 7(ri) = Ar(r2) accounting for a strength parameter A and a frequency distribution function r(r2). Note that 
T{il) refers to the effective frequency distribution of the CBR photons which interact with the system of particles 
at energy around hw. We next consider that the system of particles and CBR interacting almost resonantly with 
Lorentzian spectrum 

TT T^(S2 — (jj)^ + 1 

In view of Eq. ( ^) it follows from the Fourier transform of Eq. (|27b|) that 

dB,{t)dBj{t') = — e*'^(*-*') e-(*-*')/^- dt, (29) 
27r 

where the correlation time defines the memory time over which the stochastic function changes appreciably. From 
Eq. ( p9| ) we conclude that when considering extremely short, i.e., much less than all other times of interest (evolution 
of the part icle s ystem ) so that in a good approximation dBi{t)dBj(t') ~ 5{t — t )dt, the system is Markovian. Through 



Eqs. (27a) and (27b) the physical stochastic differential equation reads 



d\^N+CBR) = <^ -^HN+cBRdt+ U^Y, [A{n)4 + A\n)sLk) -dBin) 
[ fe=l 



dnT{n) 



■ N 

E 

.k=l 



(A{n)a.l+A\n)a.u 



dt ) \^n+cbr). (30) 



It must be emphasized that Eq. (|30|) describes the evolution of the state vector of system of N particles and CBR 
differently from the stochastic differential equation in the CSL model. The Hamiltonian H^+cbr in this equation 
describes the free evolution of the system of particles and CBR, while the two remaining terms account for the 
stochastic interaction between the CBR and partices. 

By defining both, the Wiener process dB and the operator Z depending on the CBR frequency space, the positional 
space will not be anymore privileged with respect to the momentum space, as occurs in the CSL model. We now 
proceed to the separation of the CM motion of the modified dynamical equation (BO) . The substitution of the operators 
a^,afc as position and momentum operators Pfe,qfe, permits us to express Eq. (pOl) in terms of the CM coordinates 
Q =^ Efc qfc and P = Efe Pfe as 

d\^N+CBR) = ^-'-^HN+cBR.dt + J dCt {A{n)X^ + A\n)x) ■ dB{n) 

- ^ J dnr{n) {A{n)x^ + A^n)xYdtj \^n+cbr). (3i) 

where the operator X accounting for the macroscopic object reads 

^ (iVmwQ + iP), (32) 



y2hrrwo 



while X^ is its hermitian conjugate. These operators satisfy the commutation relation [X^, ATJ] = N6ijl. As mentioned 
earlier the coupling constant of the interaction between the CBR and the system of particles defines a characteristic 
time A~^ which is associated to the net effect of the random pseudo- "potential" described by the last two terms on 
the right-hand side of Eq. (^. 

As the stochastic operator in Eq. ( ^ ) automatically acts only on the joint wave vector of the CM degree of freedom 
and the CBR cm+cbr) , the separately normalized state vectors \^ cm+cbr) and | the latter for the internal 
degrees of freedom, satisfy the equations 



CM+CBR) 



-'-HcM+cBRdt + Jdn (A(r!)xt + A\n)x) ■ dB{n) 
J dnr{n) {A{n)x'' + A\n)xf dt 



\'^cm+cbr), (33a) 
d\(f>{ri}) = -*^^(,,}|0{r,})rfi- (33b) 
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, Pn+cbr 



It sh ould be noted that the above Eqs. ( p3a| ) and ( pSq ), differently from those in the CSL model (Eqs. (17a) 
and ( |l7b| )), involve also the CBR degrees of freedom. As will be shown later, the present approach in the low 
temperature limit allows obtaining separately the normalized wave functions for the system of particles, \'^cm) and 
|0{ri})i satisfying equations similar to those in the CSL. Next, from Eqs. (||) and (^) the statistical operator 

PN+CBR = \'^n+cbr){'^n+cbr.\ rcads 

^^^j^ = -1[Hn+cbr, pn+cbr] dnvijl) [{(A(f^)Xt + At(f^)x) 

- 2 {A{n)x^ + A\^i)x) ■ PN+CBR {A{n)x^ + A^{n)x)] , (34) 

which is a prc-mastcr equation in that it contains operators from both, the N particles system and the CBR, allowing 
to calculate correlations between operators of system of particles and CBR. However, since we only have at our 
disposal the statistical properties of the CBR field, the obvious procedure is to trace over the CBR degrees of freedom, 
considered thermalized at temperature T, which leads to the reduced density operator of the system of particles only, 
containing the average number of photons of the CBR as a parameter. 

Back to Eqs. (p9|), when the following assumptions are met: i) a short correlation time Tc (<C A^^), leading to the 
Markovian approximation; ii) the interaction between the system of particles and CBR is sufficiently small (exactly the 
purpose at hand), the density operator of the global system can be written as PN+CBR{t) = PN{t)®PcBR{t)-\-Pcorrei{t), 
where the correlation term Pcorrei can be neglected |p^ . By considering the thermalized CBR density operator 
PCBR = exp (^—(3Hcbr[-^\-^)) /Tr [exp (— /3ffc'B^(A'f, A))] , with /3 = fcsT, ks being the Boltzmann's constant and 
T the CBR temperature, we find the master equation for the N-particle system 

'^ = -'-[Hm,pn]-^ J d{n)r{n){[xK,XpN] + [pNXl,x] 

+ {n)n ([Xt., [X,pn]] + [X; [Xt,p^]])} , (35) 

where pN is the reduced density operator of the N-particle system and = 1/ (exp(/3?ii7) — 1) is the thermal 

averaged photon number. 

As time goes on, it is expected that the stochastic coupling induces the iV-particle system to a thermal equilibrium 
with the CBR. By evaluating the rate of energy change between the system and the CBR we shall estimate the strength 
parameter A and improve our understanding about the nature of this stochastic coupling. In order to estimate the 
energy mean- value let us consider the mean value of a generic dynamical variable V whose equation of motion is 
obtained through Eq.(p3) as 



= _^tr ([V, Hn]pn) - 1 1 d{n)r{n)Tv { [[v, xt] . x + xt . [x,v] 

+ (n)o ([[V,Xt].,X] + [[V,X].,Xt])] p^} , (36) 

By applying Eq.(p6|) to the position and momentum variables consecutively, we observe that not only the pure 
Schrodinger evolution is modified but also the results from the CSL model, such that the equations of motion become 



d(P) 1 



dt 2 



NA(P). (37b) 



These equations lead to the results (P)t = exp (— ^iVAi) (P)s and (Q)t = exp (— ^iVAi) (Q)s, where the subscript s 
indicates the pure Schrodinger evolution: (P)^ = (P)i=o and (Q)^ = (Q)o -I- {P)t=o/M t. For V = Q^, Q • P + P • Q 
and P^ successively, the equations of motion for the mean values become, respectively, 

= -L(Q . P + P . Q) _ NA{Q^) + dnm) (1 + 2{n)n) , (38a) 

"^^^"^dt^ ^^ = ^^^^ • P + P • Q>, (38b) 

= -iVA(P^) + ^^^^^ / dnrm (1 + 2{n)n) , (38c) 



dt ' ' 2 _ 

which differ from the pure Schrodinger evolution since A 7^ 
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IV. THE MASTER EQUATION AND ITO DYNAMICS 



It will be useful to remind the conventional treatment of the problem of interaction of a N-particle system with the 
reservoir (R). Under the Hamiltonian H = + Hn + V, V being the interaction between both systems, the reduced 
density operator of the N-particle system, pN{t) = Tr^, [pN{t)], evolves, up to the second order in the interaction, 
according to the generalized master equation [|lq| 



dt 



-[HN^PNit)] 



1 



Trr 



[V,e-'^o('-''^[V,pNit)pj,]]dr, 



(39) 



where Lo(-) = [Hn + Hn, ■] is the Liouvillian operator of the free Hamiltonian. The second term in Eg. (|39|) , acting 
as a source of noise for the system and also as a sink (or source) of energy, is responsible for the irreversibility of 
the process and the loss of coherence in pN{t)- As such, the Ito calculus is justified when the stochastic terms are 
introduced into the Schrodinger equation. So, the CBR is responsible for the variation of the mean energy of the 
system and the increase of entropy. As shown by Isar et al. M, choosing conveniently the interaction term V it is 
possible to obtain Eq. ( ^Sf ) (the Lindblad form)from Eq. (|3£ y 
It is worth noting that the master equation (|35|) can be written as 



dpNjt) 
dt 



--[HN,Ps{t)] 



^S[Cn]pN{t), 



(40) 



where the superoperator 5[c„] is defined as 



S[Cn]pN = Cn.pNcl 



Pn} 



(41) 



with ci = [A / dnr{n){n)n] X'' and C2 = [A J df7r(0) (1 + {n)n)] ' X. Written as in Eq. our master 

equation resembles the Lindblad form for the decay of a mode of the eletromagnetic field inside a cavity . 

In summary, we have assumed ad-hoc that the evolution of the system of particles in its way from quantum to 
classical dynamics, under the influence of the CBR, is described by an Ito stochastic equation. However, here we 
showed that the usual master equation formalism can be viewed as a subdynamics of the Ito dynamics, without any 
need to use perturbation methods as is done in the conventional derivation. 



V. THE STRENGTH PARAMETER 



Back to the equations of motion (Bg), their solutions are 



(Q') = (Q'>.e- 



M 



t 

iVA 



({Q,P}) = ({Q,P}),e~^^*-3Ifi^ 

ZXNmhuj 



-NKt 



te 



1 - 



-NAt 



NAt 



2 J 
1 

"tva 



NAt 



1 - e 



1 



7V2A2 



1 

2^ 



-NAt\ 



-JVA^^ 



-NAt\ 



(42a) 
(42b) 
(42c) 



where 



(Q'). = (Q')o + ^(({Q,P})oi+^ 



({Q,P}), = ({Q,P})o + -(PV, 



(43a) 

(43b) 
(43c) 



The effect of the CBR temperature is present in the integral I = / drir(il) (1 + 2(n)si). It is worth noting that the 
time evolution of the operators in Eqs. ( ^ ) does not show th e ad ditive property with respect to the Schrodinger 
terms as obtained in the CSL model. As a consequence, Eq. ( [42c) ) differs from the corresponding one in the CSL 
model, Eq. (20b), because instead of the diffusion, inducing a steady increase of the mean value of the kinetic energy, 
the present model exhibits, asymptotically, thermalization due to the CBR, 
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{K) = {{K),~ K,,) e-^A* (44) 

where the equihbrium kinetic energy reads Kpq = 'iThioj'^. So, w is a characteristic frequency proportional to the 
thermalized mean kinetic energy of the CM. 

As mentioned above, in the CSL model the localization of a single particle of the system is sufficient to localize the 
whole system; as a consequence, the CM energy increases linearly with the "interaction" parameter iVAi. However, 
from Eq. we conclude that the stochastic coupling accounts for a CM energy which grows or decays exponentially 
with NKt, depending on the negative or positive value for {K)s — Keq, respectively. 

In order to estimate the strength parameter A, from Eq. (^) we assume that the relaxation time follows from 
relation {{K)s - Keq) e-^'^^^ ~ K^q, so that 

Ntr \ Keq J ' 

For a system of iV « 10^^ particles initially at room temperature the equipartition energy theorem gives a mean kinetic 
energy {K)s ~ lO^ergs. The integral I accounting for the effect of the temperature of the CBR has been estimated 
in the Appendix for (3tc <C h, with ujTc < 1. The result X ~ 1 + 2(n);^, holds for both, low- and high-frequency 
regimes. So, we find for the equilibrium energy at low-frequency regime {Hlo <C fc^T, so that {n)i^ kBT/hw), 
Keq ~ ksT ~ IQ^^^ergs. At high-frequency regime {huj » kgT), the equihbrium energy obeys K^q ^ ksT. (We 
are referring to low- and high-frequency regimes since the nowadays CBR temperature, T « 3K is assumed). Taking 
Keq at low-frequency regime (in fact, due to the In function, choosing Keq in low- or high- frequency will not change 
appreciably the value of A), and the relaxation time tr of the order of the age of the Universe, about lO^^s (what 
seems to be reasonable when considering, as obtained below, such a small coupling of the system to the CBR), we get 

A«10-38s-\ (46) 

a value to be compared with the above-mentioned coupling in the CSL model ( ~ lO"'^" cm^s~^. Thus, the parameter 
A is of the order of the upper limit of the excitation rate for nuclcons estimated by Pearle and Squires by 
comparison with neutrino-induced process. As already pointed out, such a value hardly affects the dynamics of a 
microscopic particle. 



VI. WAVE-PACKETS REDUCTION RATES 



Back to Eq. (pSi), in the CM positional representation, the density matrix pn{Q,Q') evolves according to the 
differential equation 



^Pjv(Q,Q^t) 
dt 



92 



2iM V9Q2 dQ'' 



V 



(Q-QO'- " A — + — 
pjv(Q,Q',t). 



d_ 
dQ' 



-NA 
2 



^ d ^, d 
Q h Q — 



(47) 



The first term on the right-hand side comes from the commutator in Eq. (|3^), the terms multiplied by the diffusion 
constant T> = NMAlu (1 + 2{n)^) /Ah (as well as the remaining terms, which are independent of temperature) account 
for the fluctuations (or random kicks) and for the energy changes due to the stochastic coupling, respectively. 

To analyze the wave-packet reduction rates we will not consider Eq. (|4^) in detail, since the effect of the second 
term on quantum superposition will be of much greater interest . For a brief estimation of the off-diagonal matrix 
elements of Eq. (B^) will decay exponentially as 



(Q|ps(OIQ'> -c-<*(Q|p5(0)|Q' 



(48) 



where ( = ^'(AQ)^ and (AQ)^ = (Q - Q')^. It follows from Eq. (||) that the quantum coherence of a macroscopic 
system will disappear on a decoherencc time scale 

1 1 h 



I?(AQ)2 (l + 2(n)„)iVMAt^(AQ)2- 



(49) 



Analyzing Eq. ( ^ ) in terms of the CBR temperature, it is interesting to note that in the low-temperature limit 
(nowadays universe, T ~ 3K), i.e., {n)^ 0, the number of particles N plays a crucial role in the decoherence 
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process induced by the CBR. In the high-temperature hmit, i.e, {n)^^ oo (the early universe in the present model), 
we conclude that Eq. ( ^) leads from quantum to classical physics even a system composed by a small number of 
particles. This is a key result which help supporting the assumptions considered in the present model. 

Let us now estimate the decoherence time for both, a macroscopic and a microscopic object in nowadays Universe, 
i.e, T ~ 3if . In order to compare our results with that presented in literature, we consider the low-frequency regime, 
such that Eq. (^ ) reduces to 

1 

^ I?(AQ)2 " 2A^AMfcBT(AQ)2 ' ^^^^ 

By considering a system of N (~ 10^'^) hydrogen atoms with mass M « Ig and separation AQ « 1cm, quantum 
coherence would be destroyed in r^j w l{)~^'^s. Such a value turns to be significantly smaller than the one obtained 
by GRWP, XcM ~ 10~''s , Eq . ([l9| ) , and comparable with that obtained through the linear response model of 
the Caldeira and Leggett (CL) [p3[|, where, also at low-frequency regime, td/tr « h"^ /2mkBT{AQ)^ , tr being a 
relaxation time. For the above-mentioned system of N atoms, and assuming tr w 10^^^ ^ have done to obtain 
A, Eq. (|45|), we get from CL model td ~ 10~^^s. So, Eq. (^9|), and consequently Eq. (^o|), arise from a theory that, 
despite assuring the essential character of the GRWP model, gives a more realistic value for the decoherence time of 
a macroscopic object. 

As far as a microscopic object is concerned, for example a single atom, m w IQ^'^'^g on atomic scale AQ « 10~*cto, 
we observe the persistence of quantum coherence since tjj « lO^^s. Finally, we note that when considering a tiny 
Weber bar |2^|2^, AQ « 10~^^to, at cryogenic temperatures, T « 10^^ K, we also observe the persistence of quantum 
coherence from Eq. (|4^), as should be expected. 

Back to Eq. (^), when interpreting the exponential damping factor ^ by the light of the CSL model (Eqs. ( p^ ) 
and dig])), we conclude that the strength A plays the role of a microscopic frequency hitting parameter. 

VII. THE CM AND INTERNAL MOTION 

By construction we assumed that the CBR acts only on the CM coordinates of the system of particles. Such 
assumption automatically decouples the dynamics of the collective and internal motions in the master equation (|3^). 
Next, we show that even the vector state dynamics for the CM and the internal motion decouple, as in the CSL 
model. Of course, our analysis will be restricted to the low temperature limit where, as obtained in Eq. (^9|), the 
macroscopic character of the system becomes really important due to the number of particles TV. In this limit Eq. 
( ^ ) simplifies to 

% = -'-[Hn,pn] + AX • pArXt - ^ {Xt . X, PA,} . (51) 
at h 2 ^ ' 

The stochastic differential equation for the state vector of the system of particles which leads to Eq. ( pTf ) can be 
written as 

dl^-jv) = f-^HNdt + X-dW- I^Wf • Wtj l^-s), (52) 



now with the Wiener process dWi = 0, dWidWj — ASijdt. 

The assumption made in the CSL model, that the set {ri} in Eq. (O) does not contain macroscopic variables, 
implies that the state vector for the macroscopic object factorizes as 5'Ar(]cjfc}) = '0CM(Q)0int({ri}). The additional 
assumption that the CM motion is decoupled from the internal degrees of freedom means that the Hamiltonian 
must be written as a sum of two terms, Hn — Hcm + Hint §■ Under these assumptions the Ito calculus, d^ n = 
d{il^CM4>tnt) = (dipcM)(j)int + ^C M(d(t>int) + [dipc M {d4>int) , shows that the wave functions V'cm(Q) and (/ii„t(ri). 



similarly to Eqs. ( |l7a| ) and (|l7bD , satisfy equations 

d\i;cM) = (^-'-HcMdt + X • dW - • Wdt^ \^cm), (53a) 



1 

d\4>int) = --^Hint\(f>int)dt. (53b) 
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The stochastic terms do not affect the internal structure of the system of particles, i.e., nothing changes in the 
Schrodinger dynamics of microscopic particles. It is worth noting that in the CSL model the additional assumption 
of a large enough localization width parameter (besides of an internal wave function independent of macroscopic 
variables) is necessary to decouple the dynamics of ipcM from (f>int- In fact, as shown in Ref. []lT| , a width parameter 
of order of atomic size leads to the breakdown of the translational symmetry of the system and the interaction between 
the CM and the relative coordinates {i.e., H — Hcm + Hint + V), has to be taken into account. However, in the 
present model, since we have assumed that the CBR acts only on the CM coordinates of the system of particles, no 
additional conjectures was requested about the random operator Z{^l), Eq. (p5|), to achieve the remarkable result of 
the CM decoupling from the internal motion, as if the stochastic terms in Eq. (35) were absent. The operator Z(il) 
has thus the advantage of not needing additional conjectures about the width parameter of the localization process. 



VIII. DECOHERENCE AND ENTROPY 



The decoherence process resulting from the interaction of the state vector for a macroscopic object with the CBR 
can be quantified by the rate of increase of either the linear or the statistical entropy. In terms of the density matrix, 
the statistical entropy, a measure of our ignorance, is defined as |^ Sg = — Tr (pin (/?)) (the subscript s refers to 
statistical). This definition does not require that the system be in a thermal equilibrium state. Alternatively, a good 
measure of the loss of purity for states of an evolving open system is based on the increase of the linear entropy 
(subscript I) 

Si - Tr {p - p2) _ (54) 

Next, we estimate the rate of increase of the linear entropy through the evolution of the density matrix given in 
the operator form by Eq. ([l7|). Considering a weak strength parameter (A « 0) and the state vector remaining 
approximately pure (Trp^ « Ij , up to first order in A we obtain 



5.^4P^((AQ)V^((AP^^^ 



(55) 



where ((AQ) ) and ((AP) ), obtained from Eqs. (U2a|) - ([43c|), stand for the variances of the position and momentum 
operators and can be rewritten as function of their initial values (Q)o and (P)o- 

In order to better understand the rate of increase of the linear entropy in Eq. (p5|), it is worth to compare it with 
that obtained by Zurek who used the linear response model of Caldeira and Leggett |^ (in the high temperature 
limit). With the above approximations Zurek obtained Si — 4I?((AQ)^) (for a single oscillator), so that the rate of 
increase of linear entropy (in quantum Brownian motion) is proportional to the dispersion in position coordinate only 
- the preferred observable singled out by the interaction Hamiltonian. In our approach, from Eq. ( p5[ ) we observe 
that no preferred observable emerge from the dynamic equation (^5|) (the dispersion in momentum is also present), 
contrarily even to the CSL model where the position representation is taken from the outset as privileged. However, 
for a large number of particles {N ^1), Eq. ( p5| ) indicates that the dispersion in momentum is considerably smaller 
when compared with that in position which, in this situation, emerges as the preferred observable. 

In the weak-coupling limit we integrate Eq. (|5^) replacing the general evolution in Eq. ( |35| ) by the free von 
Neumann equation to obtain 

Si^AV 
1 



3M 



l{{AQ)\ + ^^((AP)^)o W+ -i-(A{Q,P})oi^ 
r((AP)')ot^ 



(56) 



with (A{Q,P}) = ({Q,P}) — 2(Q)(P). The dispersions appearing in the Eq.(p6|) are computed for the pure initial 
state. 

Back to the preferred basis problem, we remind that Zurek considered the free Heisenberg equations for the oscillator 
operators {P,Q) and obtained the hnear entropy 21? ^((AQ)^)o + ^^pf^iuj)'^ ((A-P)^)o^ {N — 1), averaged over one 
oscillator period. So, this result corresponds only to the coefficient for the linear time-dependence in Eq.(p6|), where 
additional terms as square and cubic time-dependent behavior also take place. Such a behavior indicates that, in 
spite of the large number of particles, for large times the momentum plays an important role in the problem of the 
preferred observable because we have considered the free motion of a iV-particle system instead of a single harmonic 
oscillator. 
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IX. SUMMARY AND CONCLUSIONS 



In the GRWP model of continuous dynamical reduction of the state vector it is assumed that each microscopic 
constituent of a system of N particles is subject to a sudden collapse due to a spontaneous random hitting process 
consisting in a localization of the wave function of the particle within an appropriate range |^ . In what turns to be 
a remarkable result the localization of a single constituent of the system of particles is sufficient to localize the whole 
system. Such a spontaneous localization, considered as a fundamental physical process, induces a steady increase 
of the mean energy value of the physical system and so the increase in temperature per unit time of the universe. 
When taking into account that the age of the universe is about lO^^s, the GRWP model leads to a total temperature 
increase from the beginning of the universe of 10"'^ i^, a value claimed to be comparable with the cosmic background 
radiation (GBR) of 3K. 

In the present model for continuous dynamical reduction, also based in a stochastic differential equation describing 
a Markovian evolution of state vectors, the random hitting process in GRWP model is substituted by the intervention 
of the GBR. Such a strategy is intended to maintain (i) the principle of conservation of energy, and (ii) the claim 
that the Universe originated from the Big Bang leaving the GBR as a signature. In (i) the increase or decrease of 
the GM mean energy of the system of N particles subject to a stochastic interaction with the GBR, which acts as a 
reservoir. In (ii), taking the opposite direction to the GRWP argument (which claim that the present temperature of 
the universe comes from the increase of the total energy arising from the random hitting process), we propose that 
the GBR temperature plays an important role in the reduction of the iV-particle wave packet. So, we assumed, in 
agreement with the standard cosmology, that the Universe has originated from a hot state and cooling during its 
expansion, with decreasing mean photon energy. The Planck law for the thermal average boson number in GBR, 
indeed the best blackbody known, has recently been tested by the GOBE satellite j2^. The temperature of the GBR, 
decreasing as the mean photon energy decrease due to the cosmic expansion makes the mass of the system increasingly 
more important for the transition from quantum to classical description. On this basis one can argue that the quantum 
nature of the Universe becomes increasingly important as it is cooling. In fact, for the early Universe, the number of 
particles does not play a fundamental role in estimating the decoherence time, where higher temperatures (by itself) 
turn the system from micro to macro dynamics. However, as the Universe becomes cooler the number of particles 
becomes increasingly important. 

Moreover, the present model leads to realistic results for decoherence times. While in the GRWP model the value 
10~^s obtained for a system of particles to go from micro to macro dynamics seems to be too large, the value 10~^'^s 
here obtained for a system of N atoms in the low-frequency regime is comparable to the decoherence time obtained 
from the Galdeira-Leggett model. 

As mentioned, whereas the GRWP model requires a privileged positional space, in the present model, by construction, 
the stochastic operator acts on the GBR spectrum, carrying the same status for both, the position and the momentum 
space. The GRWP's result - the wave function collapse of a single particles induces the collapse of the wave function 
of the whole system - was obtained exactly from the choice of the position as a preferred basis. The same result 
follows from our model without the choice of the position as a preferred basis. However, It has to be mentioned that 
in spite of attributing the same status for the position and the momentum space, when analyzing the entropy under 
the process of decoherence, the position coordinate still emerges as a preferred basis when considering a system with 
a large number of particles N . So, the preferred basis is directly related to the number of particles in the system. 

Another interesting feature is that we do not claim for an additional assumption to decouple collective from internal 
motion as the required large width parameter a^^^'^ ~ lO^^cm in the GRWP model. The random operator Z{Q) here 
assumed, besides being a more conventional choice since it is associated to a reservoir (GBR), leads to the advantage 
of decoupling the GM and internal motion without additional assumption beyond that usually assumed for a reservoir. 

The random operator describing the interaction between the system and the GBR carries only one parameter, the 
strength A, instead of the two free parameters, as required in the GRWP model (a~^/^ and the mean frequency A). In 
our model, the coupling of the GBR to the system, proportional to A, corresponds to the random pseudo- "potential" 
dh of the GRWP model. As well as the parameter A in GRWP model, our A is weak enough in the sense that it 
does not affect the dynamics of a unique particle, even in the case in which its wave function is spatially spread 

Finally, we point out that the Ito equation is not derived from a physical picture of the background and associated 
scattering processes of the GBR by the system of particles. Instead of considering a particular interaction and choose 
some specific particle property sensible to the electric and magnetic field of the GBR, we approached the problem 
by modeling the interaction by a stochastic coupling, such that the dynamics could be described by an Ito equation. 
We have considered an effective strenght parameter A accounting for all kind of light-particle scattering processes. 
We also stress that our pre-master equation (34) (with respect to the particles) has still information on both, the 
system of particles and the GBR, since it contains operators of both subsystems. This approach is different from the 
usual one where for getting a master equation it is necessary to trace over the environment degrees of freedom, as is 
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done in the theories of Joos and Zeh and Caldeira-Leggett or even in quantum optics. In our model it is possible to 
calculate correlations between observables of both subsystems. However, we have get rid of CBR degrees of freedom, 
Eq. (35), just because the available information on the CBR subsystem is sparse, consisting of the blackbody radiation 
distribution function at 3K. Thus the master equation (35) expressed in the CM positional representation, Eq. (49), 
incorporates the similar equations obtained in both theories, Joos and Zeh and Caldeira-Leggett. The main difference 
between the three approaches stem in the nature of the diffusion constant (DC): In Joos and Zeh the DC originates 
from the scattering of electromagnetic waves by small objects; in Caldeira-Leggett it comes from the fluctuations 
arising from energy dissipation of the system of interest to a thermal reservoir. In our model the DC originates from 
the stochastic interaction between N particles of mass m and the CBR at temperature T. 
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APPENDIX A: CALCULATION OF INTEGRAL I 

Due to the normalized Lorentzian spectrum (Eq. (|2^)), the integral accounting for the temperature of the CBR 
readsX = 1 + 2 J dnr{n){n)u. Now, since the Planck's distribution diverges when goes to zero, the same occurs 
to the remaining integral J dilT{D,){n)Q. However, as usual, we assume that the spectrum T{fl) has its maximum far 
away from zero in order to cancel the divergence coming from {n)^. In what follows we are going to estimate under 
which conditions this approximation is valid. 

After the transformations flTc = x and 7 — h/kBTTc, the remaining integral reads 

1 /■+°° 1 1 

[X - {lot + i)] [X - (cJTc - i)] 07^ -1 ' ^^^^ 

which can be solved in the complex space through Jordan's lema, leading to the result 



2i ■ 



■+^'^-l^i%V 2^"'")br + zfl-^)lbT.-zfl + ^)l • ^^^^ 



^ 2i e7(-e+.) _i 2 "'7 + z(i - ^^)][LOT, - z(i + ^ 

It can be shown that the imaginary term coming from the above result is zero. Now, denoting 7 = p/^, where the 



parameter p is equal to ftw/fc^T whereas ^ = ojTc, the real term coming from (A2), reads 



cos(C/p)e«-l __ n 

1 £2 h 4. 2 „ (27rnn/f)21 + (47rnp2/£)2 ' ^ ' 



e«[e« -2cos(C/p)] -1-1 i^i\^+P^~ {'^T^np/£,Y] + {Airnp^/S,) 

For large n the second term of ([A^) reduces to 



p 



00 ^ 



The analysis of the above result will be restricted to the condition ^/p <C 1, with £ ^ under which the sum in 



(A4) can be disregarded (since even ^ /p <ti 1), and the first term in (A3) gives us le"''^'""^^ —1), in a way that 



the Lorentzian distribution r{Q) acts practically as a delta function {S{ri — uj)). In fact, the limit C ^ 1; leads to 
the condition uj < t^^, so that the frequency can be taken far away from zero since, as above-discussed, we are 
considering an extremely short correlation time (Markovian approximation). Under such condition it is expected that 
the lorentzian function F acts indeed as a delta function, what means that the action of the reservoir over the system 
of particles is restricted to the oscillators whose frequencies is closely related to cu. So, the problem of how far ui has 
to be from zero, in order to eliminate the divergence coming from Planck's distribution when uj 0, depends exactly 
on the lorentzian height in its maximum. Moreover, the condition £/p <C 1, with £ 1, holds for both, the low- 
and high-frequency regimes. When £ ~ 1 (so that uj ^ t^^), we get the high-frequency regime huj ^ ksT, whereas 
for £ ^ 1 even the low-frequency regime is allowed. For the latter case we have to assure that <C a; <C t^^, not 
only to get rid of the divergence arising from (n)^ , but also to hold the assumption of highly excited oscillations of 
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the CBR leading to the Markovian approximation. Summarizing, under the conditions established above we get the 
result 2^1 + 2{n)^, which holds for high- and low-frequency regime. 
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